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ABSTRACT 

An analysis of ablation processes is made for a fall-back disk with inner and outer 
radii external to the neutron-star light cylinder. The calculated ablation rate leads, 
with certain other assumptions, to a simple expression relating the inner radius and 
mean mass per unit area of any long-lived fall-back disk. Expressions for the torque 
components generated by interaction with the pulsar wind are obtained. It is not 
impossible that these could be responsible for small observable variations in pulse 
shape and spin-down rate but they are unlikely to be the source of the periodic changes 
seen in several pulsars. 
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1 INTRODUCTION 

A simple model for neutron-star disk formation from super- 
nova fall-back was described nearly twenty years ago (Michel 
1988). More recently, there has been a revival of interest in 
the subject in connection with topics such as planetary for- 
mation (Lin, Woosley & Bodenheimer 1991; Wolszczan & 
Frail 1992; Miller & Hamilton 2001), pulsar spin evolution 
(Menou, Perna & Hernquist 2001a; Blackman & Perna 2004 
), the anomalous X-ray pulsars (Chatterjee, Hernquist & 
Narayan 2000; Marsden et al 2001; Ek§i & Alpar 2005; Ek§i, 
Hernquist & Narayan 2005), and a common framework for 
the anomalous X-ray pulsars, soft gamma-ray repeaters and 
dim isolated thermal neutron stars (Alpar 2001). 

Fall-back disks in isolated neutron stars may differ from 
the accretion disks (a-disks) of binary systems in a number 
of respects. It is possible for them to be at radii beyond the 
light cylinder radius Rlc, with temperatures observable in 
the infra-red resulting from heating by X-rays and by the 
pulsar wind rather than mass transfer and viscous evolu- 
tion. A significant fraction of the mass may be in the form 
of dust grains. Computations of the optical and infra-red 
emission expected from fall-back disks have been made by 
a number of authors (Foster & Fischer 1996; Perna, Hern- 
quist & Narayan 2000), and there have also been several 
experimental searches for the infra-red excess characteristic 
of dust grains (see Bryden et al 2006). In this way, upper 
limits of the order of Mq have been obtained for disk 
meisses in a number of radio pulsars (Phillips & Chandler 
1994; Lohmer, Wolszczan & Wielebinski 2004). 

The recent observation of an excess in the case of the 
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anomalous X-ray pulsar 4U 0142-|'61 (Wang, Chakrabarty 
& Kaplan 2006) has been interpreted by these authors as 
evidence for an X-ray heated fall-back disk. The estimated 
mass of the disk is quite small, Md ~ lO~^Af0, and a fac- 
tor favouring its observation is the high X-ray luminosity 
(~ 10^® erg s~^) of the central neutron star relative to the 
typical radio pulsar. Its estimated present inner radius is 
well outside the light cylinder, ~ A.IRlc- An alternative 
analysis by Ertan et al (2007) fits the complete infra-red and 
optical spectrum of 4U 0142-1-61 by the emission of a viscous 
gaseous disk with Vi < Rlc and an extinction to that object 
identical with the estimate made independently by Durant 
& van Kerkwijk (2006a). This type of disk is much more 
complex than the one considered in the present paper and 
should properly be regarded as forming a boundary con- 
dition for the magnetosphere (see, for example, Cheng & 
Ruderman, 1991). The ablation processes considered in this 
paper would remain relevant at r > Rlc for this type of 
disk but their use to define the inner radius Vi , as in Section 
2.4, would not be possible. 

Details of the fall-back disk formation are not described 
here. Initially, the disk must have been internally ionized and 
viscous, with mass and angular momentum transfer rates 
intrinsic to the formation process. We assume that the es- 
sential ideas of the evolutionary path described by Menou, 
Perna & Hernquist (2001b) are valid. The thermal ioniza- 
tion instability, in which free electron recombination causes 
a sudden decrease in opacity, occurs first at the outer disk 
radius but then propagates very rapidly inward. We accept 
the conclusion of Menou et al that a neutral, passive, gaseous 
disk is then formed, with no obvious source of internal vis- 
cosity. The power input to the disk, considered here in Sec- 
tions 2.1 and 2.4, is not adequate to maintain temperatures 
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above the thermal ionization instabihty temperature and so 
does not materially affect the onset and propagation of in- 
stability except, possibly, to delay it. 

To summarize, the present paper assumes that the fall- 
back disks with which it is concerned arc passive, thin and 
have internal temperatures well below the ionization insta- 
bility temperature. Hence, following Menou et al, there is 
negligible internal ionization and the viscosity of a-disks is 
not present. We shall be concerned only with the interac- 
tion of the pulsar wind with the disk surface regions so that 
it will be possible to neglect the complex evolutionary dy- 
namics of neutral dust and debris disks (see, for example, 
Frank, King & Raino 2002). There is certainly interaction 
between neutral dust grains but to a good approximation 
we can treat the disk as a set of annuli each having a ra- 
dius r and Kepler velocity vk- Waxping of the disk then 
occurs, in general, if torque components are r-dependent. It 
is assumed here that, in the general case, the disk plane at 
formation has a finite uniform tilt angle l3, but that the disk 
angular momentum is more nearly parallel rather than an- 
tiparallel with the neutron-star spin. We consider only disks 
that are external to the light cylinder, n > Rlc, and regard 
their effect on the pulsar magnctosplicrc as a partial termi- 
nation of the pulsar wind rather than a boundary condition 
to be imposed on the solution for the magnetospheric fields. 
This limitation on n means that the equivalent force on 
the disk derived from the Lense-Thirring effect (see, for ex- 
ample, Wilkins 1972) is negligible in comparison with those 
found in Section 3.1. The Robertson-Poynting effect can also 
be neglected. The reason is that the equivalent force on the 
disk depends on its Kepler velocity, vk >C c, but the torque 
components derived in Section 3.1 from the azimuthal part 
of the pulsar wind arc independent of this factor. 

The present paper is addressed principally to three 
problems. These are: (i) an analysis of disk ablation pro- 
cesses and estimates of the ablation rate and of the infra-red 
luminosity; (ii) calculation of the alignment and precessional 
torques acting between disk and neutron star; (iii) the ques- 
tion of whether or not disk alignment, counter-alignment or 
precession could give rise to time-variable phenomena that 
are observable over intervals of no more than several years. 

The last of these has been prompted by the observation 
of periodic changes iu timing residuals and in avcragc-pulse 
characteristics (Stairs, Lyne & Shemar 2000; Shabanova, 
Lyne & Urama 2001; Haberl et al 2006) of several pulsars 
or isolated neutron stars that appear to be most simply ex- 
plained by Eulerian precession (for a recent review of this 
interpretation, we refer to Link 2006). If the existence of 
Eulerian precession were established unambiguously, there 
would be far-reaching consequences for our present under- 
standing of neutron-star internal structure. For this reason, 
we attempt in Section 4 to see if residues of fall-back disks 
could exist in these cases with inner radii small enough to 
give variability on the observed time-scales. 



2 DISK IONIZATION AND ABLATION 

2.1 Momentum densities in the vacuum solution 

In order to study the interaction between a disk and the pul- 
sar wind we shall use inertial frame cartesian coordinates in 



which the neutron-star spin angular velocity f2 is parallel 
with the 2;-a:xis. Where convenient, this frame is also repre- 
sented by spherical polar coordinates r, 9, cj). The magnetic 
dipole moment is /it = BoR^ , where R is the neutron star 
radius and Bo the equatorial surface field. It is at an angle 
^ with $7. There has been much computational work on the 
aligned neutron-star magnetosphere, but in the general case 
with ^ 7^ 0, the work of Spitkovsky (2006), based on relar 
tivistic force- free electrodynamics, seems to be the only pub- 
lished solution. Thus our present investigation is based on 
the Deutsch vacuum solution for the electromagnetic fields 
of a rotating neutron star in the form given in recent papers 
by Michel & Li (1999) and Eksi & Alpar (2005). We note 
that there appears to be some disagreement about a small 
number of near-field terms in this solution (see also Ferrari 
& Trussoni 1973; Good & Ng 1985) but these are not of pri- 
mary concern for the torque calculations made here. Michel 
& Li also observe that the static radial electric field term, 
which remains significant at the light cylinder, is in prin- 
ciple undetermined because the total electric charge of the 
star and magnetosphere is an unknown quantity dependent 
on the details of formation. Prom the electric displacement 
D and magnetic flux density B given by this solution in the 
inertial frame at r > Rlc it is simple to write down the 
spherical polar components of the momentum density, 

P=T^DxB. (1) 

47rc 

For the torque calculation described in Section 3.1, we re- 
quire them averaged over the pulsar rotation period. They 
are: 

, , M^n^sin^^ / 2 a\ /n\ 



(Pe> = 0, 



(P<^>> = 



H^fl cos^ ^sin 9 



fi 



sin ^sin^. 



(3) 



(4) 



The first term in equation (4) is derived from the product 
of two static irrotational field components and hence gives a 
solenoidal contribution to (p,^). Taken at face value, this rep- 
resents a purely azimuthal component of momentum density 
whose presence in the near field of an aligned rotating mag- 
netic dipole appears not implausible. Integration to obtain 
the rate of outward transfer of angular momentum across a 
sphere of large radius r gives the vacuum torque. 



/I P27V 
d{cose) / d0(r X p), 

of which the ^-component is the spin-down torque. 



■sin 



(5) 



(6) 



3c3 

with (Tyx) = (Tyy) = 0, as expected in the inertial frame. 
Naturally, identical results are obtained by integrating the 
time-averaged torque, formed directly from the Maxwell ten- 
sor, over the same surface. (In a frame of reference corotating 
with the star, the time-averaged components of Ty are all 
finite and r„ oc (f2 x /x) x /x; see Davis & Goldstein 1970. 
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Thus the direction of CI changes in this frame. Therefore, 
in general, it also changes relative to the dipole moment ji 
which may move in this frame during pulsar evolution in an 
uncertain manner.) 

The relation between wind momentum density and 
torque expressed by equation (5) is no more than classical 
mechanics and remains valid, with {Ty) replaced by the true 
torque (F), whatever the degrees of freedom contributing to 
p^. Thus (p^) is necessarily finite in a physical pulsar wind, 
in which the electromagnetic fields arc loaded by a density of 
relativistic charged particles. However, it is unclear whether 
or not (pe) would remain zero in a real system: a finite value 
would not contribute to (F), whose x and y components 
would necessarily be zero in the incrtial frame. The effect 
of these terms on the interaction between disk and pulsar 
wind does not seem to have been considered previously and 
is significant in the case of thin dust and debris disks similar 
to those found in the major planets or in prc-main-sequence 
stars (see, for example, Beckwith et al 1990) which present 
a very small profile to pr. 



2.2 Composition of the pulsar wind 

A pulsar luminosity, observed outside the light cylinder, can 
be broadly divided into a wind component and black- 
body radiation from the neutron-star surface in the form of 
an X-ray component jCbol- The wind luminosity, which can 
be estimated from Q and the observed spin-down rate f2, is 
further divided into a Poynting flux of electromagnetic fields 
and a relativistic particle component, C^i = £-em + -Cp. In 
previous work on disk ablation by the pulsar wind. Miller & 
Hamilton (2001) assumed Lp to be the larger wind compo- 
nent, and of baryonic composition. 

However, this view of wind composition appears to be 
at variance with theories of particle acceleration at r < Rlc 
that are widely accepted to bo at least qualitatively valid. 
In the corotating frame of reference, there may exist finite 
electric-field gaps in the open magnetic flux-line regions of 
the magnctosphere. The total current flow through these 
is limited by the Goldreich- Julian charge density, acj = 
—fl ■ B/27rc, where B is the magnetic flux density. Thus 
the total rate of baryon loss estimated from the polar-cap 
area wR^fl/c intersected by open magnetic flux lines and 
from the polar-cap Goldreich-Julian charge density can be 
at most of the order of 



(7) 



where /i is the neutron-star magnetic dipole moment. But 
the size of the finite-field gaps is limited by intense electron- 
positron pair formation with the consequence that the 
Goldreich-Julian current density acjc must contain both 
ion and positron components. In polar-cap models, particles 
are typically accelerated to energies of the order of 10^ GcV 
per unit charge, though greater energies may be reached 
in outer-gap models. It is also widely accepted that, as a 
consequence of the gap-limitation process, large numbers of 
electron-positron pairs are formed in open magnetic flux- 
line regions external to the gaps. Thus the total rate of pair 
formation can be given as TZep = nTZi, but even its order of 
magnitude, k ~ 10^~^, is not well known. The energies of 
these pairs are several orders of magnitude lower than the 



typical gap energy of 10"^ GeV. It is obvious that the particle 
luminosity, estimated in the region of the light cylinder from 
equation (7) and from the gap energy, must be jCp ^ jCw in 
most cases. For example, from the ATNF Pulsar Catalogue 
parameters given for PSR 1828-11 (Manchester et al 2005), 
we find 7^6 = 8x 10^^ s"^ giving an estimated jCp « 1.3x 10^^ 
erg s~^, to be compared with £„, = 3.6 x 10^'' erg s~^. This 
is consistent with the usual assumption (see, for example, 
Melatos & Melrose 1996) that the Poynting flux luminosity 
£em is the larger component of the wind at the light cylin- 
der. We shall assume that this is maintained at disk radii. 



2.3 Ablation processes 

In this paper, we assume that the presence of a partially 
ionized disk at r > Rlc, well outside any possible Alfven ra- 
dius, has little effect on flelds and particle fluxes at r < Rlc 
and can be seen as a partial wind termination rather than 
a modified boundary condition. The sequence of processes 
in disk ablation is described in this Section. We show first 
that the thermal X-ray flux Cboi ionizes the outer regions of 
the disk surface and so makes possible the conversion of the 
Poynting flux £em to ion or proton kinetic energy. Ablation 
is then the result of several possible secondary processes, 
including neutron production in the interaction of the accel- 
erated ions and protons with disk nuclei. 

It is necessary to find the conditions under which a disk 
consisting of neutral hydrogen, atoms of mass mA,z, and 
dust grains of radius a has significant ionization. We shall 
assume that the disk has small thickness h and that it is 
internally neutral, with a temperature well below the ion- 
ization instability temperatures found by Menou, Perna & 
Hcrnquist (2001b) for metal-rich compositions. But the sur- 
face of a thin disk must have a finite temperature owing to 
its interaction with the particle component of the wind mo- 
mentum density p^. As a result, a diffuse low-density surface 
region is present which is exposed to the X-ray component 
^boi, principally photons of average energy £ w 2A;sTs de- 
termined by the neutron-star surface temperature Tg. The 
gravitational restoring force on an atom, acting toward the 
plane of a thin disk, in terms of the Kepler angular frequency 
Qk, is mA,zClj(W at height w above the plane of the disk. 
Thus the depth of the diffuse surface must be of the order 
of 
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(8) 



typically ~ 10^ cm for hydrogen at T = 10^ K and for 

Q.K = 10"^ rad s-^. The mean life of a neutral hydrogen 
atom in this region against photoelectric ionization is 

..e=^^, (9) 

where A is the mass attenuation length for low-energy pho- 
tons in hydrogen and rriH is the hydrogen atom mass. For a 
typical radio pulsar, £ = 10^ eV equivalent to Cboi = 8 x 10^^ 
erg s~^ for a neutron star radius -R = 10^ cm. The Particle 
Data Group (Yao et al 2006) give A = 1.0 x lO"'* g cm"^ 
so that even at a radius r = 10^^ cm, the mean life is very 
short, Tpe = 15 s. For higher atomic numbers Z, lifetimes 
against ionization are of the same order. 

The extent of dust grain sublimation is uncertain. Dust 
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grains exposed to £toi have an estimated equilibrium tem- 
perature Tboi = Tsy^R/2r « 1300 K, provided details of 
grain composition, emissivity and shape are neglected. Sub- 
limation rates may be significant at this temperature which 
is also approximately equal to the K-band Wion's displace- 
ment law temperature. But its dependence on Tg and r make 
it impossible to give any general statement. This estimate 
assumes a grain radius a ^ 10~* cm such that almost all 
incident photons interact. The mass stopping power for low 
energy electrons has boon tabulated by Seltzer & Borgor 
(1982). It is so large that there is no doubt that the en- 
ergy of the emitted photo-electrons is contained within dust 
grains of these radii. Smaller grains, with incomplete reten- 
tion of photo-electron energy, have lower radiative equilib- 
rium temperatures. On the basis of the Dulong and Petit 
law, the thermal equilibrium temperature is reached within 
times short compared with the Kepler period of the disk. 
The ionization time Tpe is also smaller than the Kepler pe- 
riod leading to the conclusion that, for interesting intervals 
of jCboi and r, matter in the very diffuse surface regions of a 
disk has a high degree of ionization independent of the state 
of the dust component. Neutron star blackbody photons 
therefore have an important effect on the state of the disk 
but their energies, £ ~ 10^ eV, are several orders of mag- 
nitude too small to eject protons (a proton with 10^ eV/c 
momentum has negligible kinetic energy; see also Miller & 
Hamilton 2001). 

Following the discussion of the particle component of 
£to in Section 2.2, we shall assume that the fields at the 
disk radius, though not necessarily at radii that are many 
orders of magnitude greater, satisfy the ideal magnetohydro- 
dyriamic condition cE = — vo x B, where vo is the particle 
velocity (see, for example, Melatos & Melrose 1996; for a 
review of pulsar wind nebulae, see Gaensler & Slane 2006). 
To the extent that this is true, the further conditions E < B 
and E ■ B = are also satisfied. In this case, the qualitative 
details of particle acceleration as a result of the interaction 
of the Poynting flux in the wind with the diffuse ionized disk 
surface can be analyzed easily by making a Lorentz trans- 
formation of velocity 



from the inertial frame of Section 2.1 to a frame with fields 
E', B' in which E' = (see, for example. Landau & Lifshitz 
1962). In this frame, B' is parallel with B and of magnitude 
B' = s/B'^ — E^. A proton formed by photoelectric dissoci- 
ation of neutral hydrogen has an initial velocity (neglecting 
Kepler and thermal velocities) in this frame of — v' and its 
subsequent orbit is a circle of radius w'/wb in a plane perpen- 
dicular to B', with angular frequency ujb = ev'B' /cq' where 
q' is its momentum in this frame. Thus its time-averaged 
velocity in the original inertial frame is simply v'. Although 
v' and Vo are not, in general, exactly parallel, it is broadly 
correct to say that the proton is transported forward, with 
the existing particle component of into the body of the 
disk. As this process continues, the Poynting flux is nec- 
essarily converted to kinetic energy. The energy transfer is 
linearly dependent on rest mass and thus is almost entirely 
to protons or partially ionized atoms. The penetration depth 
is very roughly defined by equating £em/47rr^c with the par- 



ticle pressure, as for the coUisionless shock present in pulsar 
wind nebulae (Rees & Gunn 1974). 

In this elementary analysis, the time-dependence of the 
E and B fields in the inertial frame has been ignored. This 
appears a reasonable approximation because, at disk radii 
r ~ Rlc, the proton lob is several orders of magnitude 
greater than f2. We emphasize that the diffuseness of the 
disk surface is the important factor in its interaction with 
the Poynting flux component of the pulsar wind. 

Disk ablation, which we define as baxyon loss, can occur 
principally through three processes. The initial assumption 
made here is that the most important is the production of 
neutrons in strong interactions between the accelerated pro- 
tons or with heavier nuclei. The nuclear interaction cross- 
section for baryons is defined primarily by the nuclear ra- 
dius. Thus the mass attenuation length is A;, w 100 g cm~^. 
A second source of loss is the possibility that, with the re- 
flection of E and B fields that is neglected here, protons and 
ions may be transported outward and away from the surface. 
Finally, the secondary protons formed in relativistic proton 
interactions inside the disk can also acquire a momentum 
component parallel or antiparallel with B which allows them 
to move out of the disk in a time small compared with the 
neutron star rotation period. These processes are all inde- 
pendent of the state of the dust component in the disk. The 
order of magnitude ablation rates obtained from the first 
of them are very different from those of Miller & Hamilton 
(2001) which, apart from a solid angle factor, were simply 
equal to jCw divided by the gravitional potential energy of a 
disk proton. 

Neutron production in strong interactions in the body of 
the disk has characteristics similar to those found by Agosteo 
et al (2005) who have measured the neutron spectra given 
by the interaction of 40 GeV/c protons with a copper target. 
(That the nuclear charge Z = 29 exceeds the disk average 
is unimportant for this purpose.) There are two groups of 
neutrons. Qualitatively, forward-directed neutrons interact 
with disk nuclei, as do secondary protons and mesons, to 
produce a hadronic cascade with length scale determined by 
Xb- The more numerous neutrons are from decay of target 
nuclei and are emitted isotropically with average energy ^ 3 
MeV. These neutrons interacting with disk nuclei are either 
moderated by elastic scattering or undergo (n, 7) capture. 
The photon component of the cascades, mostly from tt" me- 
son decay, can be a source of neutrons through excitation of 
the giant dipole state in nuclei, but this is unimportant com- 
pared with direct strong-interaction production. Provided 
the disk is optically thick, the photon energy passes to ther- 
mal degrees of freedom through Compton scattering and 
photoelectric absorption. Very qualitatively, we can assume 
that, in a high-energy cascade, protons of momentum less 
than 1 GeV/c do not interact further but come to rest after 
losing energy by ionization in accordance with the Bethe- 
Bloch formula (see Yao et al 2006) . For these considerations, 
it does not matter whether the target nuclei in the body of 
the disk, where the E and B fields are much reduced by 
conversion of the Poynting fiux to kinetic energy, are free or 
still bound as atoms in dust grains. 

For disks that are optically thick for high energy pro- 
tons, satisfying ph S> A;,, where p is the matter density, dif- 
fusion of the low-energy group of neutrons backward to the 
surface is the important mechanism of ablation. Obviously, 
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this process is dependent on the angle x between and the 
disk outward normal unit vector nx. The neutron produc- 
tion rate grows approximately exponentially with depth, the 
scale length being ~ A;, sin x, but the probability of diffusion 
to the surface without capture decreases more rapidly. It ap- 
pears that only the primary and a small number of secondary 
target nuclei in a cascade are effective sources of neutrons. 
On this basis, and the neutron production measurements of 
Agosteo et al, we shall adopt a number A/", in the interval 
1 < A/" < 10, as the mean number of neutrons escaping the 
disk per proton accelerated by the incident Poynting flux 
as described above. We shall also assume a constant value 
£p = 10 GeV for the accelerated energy. 

Neutron diffusion is more probable from disks that are 
not optically thick to high-energy protons, for example, 
those with ph ~ lOXb- Thus the ablation rates of these can 
be higher, as noted by Miller Sz Hamilton (2001). There is 
also the interesting possibility in these instances that the net 
momentum transfer normal to the disk may be antiparallel 
with to the normal component of (p^). The arguments for 
this will be given in Section 3.3 with a brief discussion of 
some consequences that might be observable. 

It must be obvious that the analysis of ablation pro- 
cesses given here is intended to be no more than qualitative 
in the extreme. In particular, there is no rigorous solution for 
the interaction of the E and B fields with the diffuse outer 
regions of the disk. One difficulty is that the value of B at 
the light cylinder varies by four orders of magnitude for the 
specific neutron stars considered here. Thus the description 
of Poynting flux conversion to particle kinetic energy given 
here is more obviously valid for the smaller values of the 
orbit radii v' /ub than for the larger values present in long- 
period neutron stars such as 4U 0142-1-61. But these failings 
can be viewed less seriously when it is remembered that the 
properties of the physical pulsar wind, even in the absence 
of a disk, are not well known. The present analysis attempts 
no more than to show that, on interaction with the disk, 
much of the Poynting flux is irreversibly converted to pro- 
ton kinetic energy and that the average proton momentum 
is roughly parallel with the Poynting vector. 



2.4 The ablation rate and disk luminosity 

We consider first the effect of the time-averaged azimuthal 
pulsar wind component (p,^) on the disk. The coordinate 
system adopted is that defined in Section 2.1. and the disk 
is assumed initially plane, with uniform tilt angle f3 and rec- 
tilinear line of nodes in the xy-plane at an angle 7 with the 
X-axis. For integration over the area of the disk it is conve- 
nient to transform from spherical polar coordinates 6, (f> to 
an azimuthal angle ^p, defined in the disk plane and with 
respect to the line of nodes. The angular relations satisfied 
by points on the disk are: 



cosx 
cos 9 

sm{<p — 7) 

cos(c^ — 7) 



— cos(0 — 7) sin 13, 

sin B sin tjj, 

cos P sin ip 

•\/l — sin^ /? sin^ ijj 
cos ip 
— sin^ /? sin^ ip 



The power input per unit area of disk is —c^{p^ ■ n±). For a 
disk annulus of radius r and width 5r, the total power input 
to both sides is. 



5W 



= rSr 

Jo 



dipc \{p4, ■ nx)| . 



(12) 



It will be convenient to eliminate n in favour of Cw = f2(r„z) 
using equations (2), (4) and (6). Given that 



(13) 



dtp sin |cos((/> — 7)! =4, 

equation (12) can be re-expressed as, 
3c^£,w sin /3 



5W = rSr 



2nQ^ sin^ ^ 



(14) 



The first term on the right-hand side of this equation derives 
from the solenoidal component of (p^) which was mentioned 
in Section 2.1. But its singularity at ^ = is clearly an 
artefact of the unphysical vacuum solution for the E and B 
fields and for disks external to the light cylinder with radii 
Ti < r < To, we need retain only the final term in equation 
(14). Integration gives the total estimated power input to 
the disk from the azimuthal wind component, 

W = ^RlcjCu, sin/3 (--—)■ (15) 
27r Kn To) 

The rate of baryon loss per unit area of disk at radius r is, 

sin/3 



B 



(16) 



(11) 



Evaluation of equations (15) and (16) for the estimated disk 
parameters of 4U 0142-1-61 (Wang et al 2006) is of interest. 
The wind luminosity derived from the ATNF parameters 
(Manchester et al 2005) is 1.2 x 10^^ erg s~^ and the inner 
and outer radii are 4.7Rlc and 16-Rlc, respectively. The 
total power input to the disk is then W = 9 x 10^° sin /3 
erg s"'^. At the inner radius, the power input per unit area 
is 5 X lO'^sin/? erg cm^^ s^^ which is large in relation 
to the gravitational potential energy per unit area of disk 
divided by the to ~ 10^ yr lifetime of the neutron star, 
GMT, /ton ~ 5 X 10® erg cm~^ s~\ where M is the neutron 
star mass and S « 1.6 x 10'* g cm~^ is the mean mass per 
unit area found from the disk mass and radii given by Wang 
et al. Thus the inner edge power input from the wind is large 
compared with any possible rate of thermal energy release 
by viscous evolution of the disk. Assuming re-radiation ex- 
clusively in the infra-red, this power input corresponds with 
a blackbody temperature of T 960 K, which is lower than 
the temperature of 1200 K found by Wang et al. This is not 
unreasonable owing to the very large value of Choi ^ C-w 
for this neutron star. We have assumed the disk thickness 
to be small and uniform. But it is possible that it has some 
radial dependence, possibly of the form h cx r^^'' which is 
characteristic of point-source illumination of an a-disk (see, 
for example, Pringle 1996), and that some of the power in- 
put to the disk is from Cboi as in the model fit of Wang et 
al. The model adopted by these authors must include some 
assumption about the absolute thickness of the disk, in ef- 
fect, the solid angle it subtends at the neutron star. This 
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does not appear to be stated, but it is also of interest that 
their fit involves a very large X-ray albedo, rid = 0.97. Fi- 
nally, we have to accept that non-disk contributions to the 
the infra-red luminosity may also be significant. It is possi- 
ble that the recent observation of short timc-scalc downward 
fluctuations in the K-band intensity of this object (Durant 
& van Kerkwijk 2006b), not readily accommodated within 
the passive disk model, may be evidence for these. 

Evaluation of equation (16) at the present value of O 
gives an estimated baryon loss rate, at the inner radius, of 
3 X 10^° cm-2 s"^ for A/" = 10 and £p = 10 GeV, wjiich is 
negligible in comparison with the mean disk density E. This 
is an essential result because the ablation rate must have 
been considerably higher at earlier times following neutron- 
star formation. It is more interesting to integrate the abla^ 
tion rate over the life of the star, assuming simple spin-down 
given by equation (6), and to compare the total baryon loss 
per unit area at the inner disk radius with E. We have. 



/' 

Jo 



Bdt : 



3c7 (MsinP 



(no - n) , 



(17) 



where I is the neutron star moment of inertia and fio its spin 
angular frequency at formation. This assumes that, at earlier 
times, the disk extended inward to smaller radii but has 
since suffered ablation to its present value of rj. It provides 
an estimate of a combination of the unknown parameters 
£p and N in terms of quantities accessible to experimental 
measurement, 



ScItuh 



(18) 



The comparison for 4U 0142-1-61 gives a value for the right- 
hand side of equation (18) of 3.4 x lO'^ erg~^ s~^ or 5 x 10** 
rad s~^ GeV~^. Insertion of the values we adopted in Sec- 
tion 2.3 for Af and £p then leads to sin /3 = 5 x 10'* which 
is obviously about two orders of magnitude too large. But 
given the considerable uncertainties both in our parameter 
values and in those estimated by Wang et al, the discrep- 
ancy is not disturbing. If the disk parameters of Wang et al 
are accepted, the comparison suggests smaller £p and larger 
Af than the values adopted in Section 2.3. Possibly, the two 
processes of baryon loss listed but not considered there may 
be more important than we assumed. We have also neglected 
any increase or decrease in disk radii caused by angular mo- 
mentum transfer from the neutron star via the torque com- 
ponent T± given by equation (24) in Section 3.1. 

Apart from the uncertainties in A/" and £p, we emphasize 
that there are at least two principal sources of uncertainty 
in our calculation of W. Our treatment of the wind-disk in- 
teraction in Section 2.3 is elementary and does not give any 
estimate of the reflection of wind energy. There is also the 
uncertainty in our assumption that pe makes no contribu- 
tion because its time-average vanishes. This may not be too 
unsatisfactory for torque estimates, as in Section 3.1, but 
may nevertheless introduce an error in calculating the total 
power input. 



3 DISK-PULSAR TORQUES 
3.1 Calculation of the torque 

Calculation of the torque components acting on the disk 
follows that for 5W given by equations (11) and (12). It 
will be convenient consider an annulus of width 5r and to 
resolve the torque acting on it in the plane of the disk into 

components STp and 5Ta that are, respectively, parallel with 
and perpendicular to the line of nodes. The components are. 



6Tp = —cr Sr 



and 



/ dV'sinV'(p^) |cosx| 
Jo 

(cosx+ |sgn(cosx)) , 



(19) 



6Ta = cr 6r / d-ip cos tp {p,j,) \cosx\ 
Jo 

(cos X + g sgn(cos x)) , (20) 

formed from the momentum flux perpendicular to the plane 
(neglecting any wind-energy albedo) and with the assump- 
tion that the disk is optically thick and re-radiates photons 
as a Lambertian surface. Use of the time-averaged azimuthal 
momentum density in these expressions is valid because the 
neutron star rotation period is several orders of magnitude 
smaller than the Kepler period of the disk. Thus it is per- 
missible to treat the annulus, though not the whole disk, as 
a rigid body. Evaluation of the first of these integrals gives 
a zero disk precession torque, SVp = 0. This cancellation of 
contributions from the four quadrants of ■!/> is a consequence 
of the vacuum- field form of (p,/>), as given by equation (4), 
and it is by no means obvious that it would hold for a phys- 
ical pulsar wind. Integration of equation (20) gives a finite 
torque which acts to align the angular momentum (5Ld of 
the disk annulus with the neutron-star spin. 



SVa = cr 5r 



in which the integral. 



^/2 



cos^ V 



-^/l — sin^ (3 s\v? tp 



(21) 



(22) 



is a slowly varying function of 0, 0.66 < la < 0.79. 

The torque component perpendicular to the plane of 
the disk is, 

5r± = cr^'Sr / d-ip{p4,) |cosx| 
Jo 

(cos l3 sin 6 + sin cos 6 sin(<i6 — 7)) . (23) 
Evaluation of the angular integral using equations (11) gives, 
3c^^ 



5T± = cr Sr 



(2/3cos/3(l-Fsin^ /?)) 
1 



877^3 sin^ ^ 
^cos^^+^^ + ^Jsin^Cj. (24) 
We shall assume that the angular momentum of a fall-back 
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disk is likely to be more nearly parallel, than anti-parallel, 
with the neutron star spin, in which case the interaction 
increases 5Ld and moves the disk outward to larger radii. 



3.2 The efFect of torques on the neutron star 

The motivation for considering this aspect of the interac- 
tion is that the angular momentum L of the neutron star 
can be several orders of magnitude smaller than that of the 
disk. In the case of 4U 0142-1-61, for example, the disk pa- 
rameters estimated by Wang et al, including its very small 
mass Md ~ W'^Mq, give Ld = 2.0 x 10'*'^ g cm^ s"^ which 
is between one and two orders of magnitude larger than L. 
Consequently, the changes in direction of both L and Ld 
can be significant. Thus the Euler equation for the neutron- 
star angular momentum in the absence of a disk, L = F, is 
replaced by, 

L Ld = (f - Td) (Fd) , (25) 

for the whole system, in which Fd = Fa -I- Fp -I- Fx is the 
integrated torque acting on the disk, and F is the modi- 
fied non-disk torque. For hypothetical disks with Vi ^ Rlc, 
the terms in equations (21) and (24) that are significant 
contributors to the torque would be those derived from the 
long-range component of in equation (4) . But the vacuum 
torque (F„) can also be obtained from this term, as could 
the true (F) if the true p,^ were known (see equation 5). This 
suggests, but does not prove, that the first, bracketed, term 
on the right-hand side of equation (25) approaches (F) in 
this limit, the disk torque being cancelled by terms in (F). 
However, there is no reason to suppose that this cancella- 
tion is exact at disk radii r ~ Rlc of physical interest where 
the near-field terms in equation (4), (21) and (24) are not 
negligible. 

The terms in Fd, including any non-zero Fp, are all 
finite in the inertial frame when averaged over the neutron 
star period. Thus the direction of L must change during the 

evolution of the system, but at a rate that is unlikely to 
be of observational interest. To confirm that this is so, we 
can consider PSR 1828-11 as an example and suppose that 
L is at a small angle u with the total angular momentum 
L -I- Ld and precesses about it with angular frequency Qp. 
The torque required is Fd ~ Lufip — 2 x 10^^ g cm^ s~'^ for 
a 500 day precession with amplitude lO^'^, several orders 
of magnitude greater than the spin-down torque which is 
£yj/n = 2.3 X 10^^ g cm^ s"^ for this pulsar. 



3.3 Alignment and precesssion of the disk 

The initial evolution of a very thin disk with small E and 
Ld <^ L can be obtained directly from equation (25), to 
the extent that we can neglect the effect of changes in its 
geometrical form on calculations of Fd. The solution is, 

Ld = Fx 

LdP = -Ta 

Ld sin p-y = Fp. (26) 

Prom a hypothetical initial state of a plane disk with 13 ^ 0, 
the alignment torque Fa reduces /? at a radius-dependent 



initial rate given by 

A ^ _sra ^ 1 /3C^\ 

^ SLd 27rS(GM)i/2r5/2 [snnj 

(4/aSin^/?-h|7rsin/3) , (27) 

in which only the long-range term in equation (21) has 
been retained. The disk mass per unit area E in equation 
(27) is the local value, not the mean S. There is initially 
no precession for the vacuum-field pulsar wind assumed in 
Section 3.1. But the r-dependent alignment rate destroys 
the planar form of the disk and therefore introduces a fur- 
ther torque contribution derived from the radial momentum 
flux c{pr) + Cboi/^Tvr'^c. This is the torque concerned in the 
Pringle instability (Pringle 1996; see also Petterson 1977, 
Prank et al 2002). Initially, it is a contribution to the pre- 
cessional component Fp, but further evolution of the shape 
and motion of the disk is more complex. It is worth adding 
that any finite torque component derived from equation (19) 
as the consequence of a momentum density (pe) in a 
physical wind, rather than the vacuum field solution used in 
Section 2.1, would also induce precession. 

In relation to disk alignment and precession, the inter- 
esting question is whether or not these processes are suffi- 
ciently rapid to produce time-varying phenomena that are 
observable within periods of no more than several years. Por 
phenomena involving the body of a 1O~^M0 disk, the an- 
swer must be in the negative because Fd is perhaps an order 
of magnitude smaller than the spin-down torque F^^ yet, as 
we have noted in Section 3.2, the disk angular momentum 
Ld is at least of the same order as L. 

But the inner edge of a disk must have small local val- 
ues of E owing to its continual ablation by baryon loss, given 
by equation (16). In principle, the alignment rate given by 
equation (27) for these radii could become observably fast. 
There is also the likelihood, mentioned briefly in the penul- 
timate paragraph of Section 2.3, that for some intervals of 
small E, the net momentum transfer normal to the disk is 
antiparallel with the normal component of (p^) . 

We can see how this arises by referring back to both the 
description of hadronic cascade development in Section 2.3 
and the torque calculations of Section 3.1. The integrands 
in equations (19) and (20) both contain the rate of momen- 
tum transfer normal to the disk and per unit area. Both 
expressions are based on two assumptions; that the disk re- 
radiates as a Lambertian surface and that it does so from 
the surface of incidence with which the wind interacts. The 
second assumption is certainly valid if the disk is optically 
thick for the photon spectrum emitted and also has suffi- 
cient thickness, at least of the order of E ~ lO'^ \b \ cos x\, 
that the hadronic cascades terminate much nearer the sur- 
face of incidence rather than the far surface. But the transfer 
of energy in a cascade from the primary particle to thermal 
degrees of freedom and to low-energy photons of perhaps 
lO'^ — lO** eV occurs predominantly in the late stages of its 
development. Thus for E ~ 10Af,| cosxl, the wind heats the 
fax surface of the disk, not the surface of incidence. The ef- 
fect of this is to change the sign of the sgn{cosx) terms in 
equations (19) and (20) and, unless x is relatively small, to 
change the sign of ST a, transforming it from an alignment to 
a counter-alignment torque. There is a second factor having 
the same effect. The production of low-energy neutrons re- 
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ferred to in Section 2.3 also occurs mainly in the late stages 
of cascade development. For these E, the baryon loss rate 
per unit area much exceeds that given by equation (16) and 
is mainly through the far surface of the disk. This is sig- 
nificant because massive low-energy particles are effective 
carriers of momentum. 

Validation of these qualitative arguments would require 
very extensive Monte Carlo calculations of cascade develop- 
ment. But although these have not been made, wo suggest 
that there is a strong likelihood that a counter-alignment 
torque acts on the inner edge of a disk. Again, we shall not 
attempt to consider the complex and detailed evolution of 
a disk under those torque components, but restrict our dis- 
cussion to the question of whether or not they could cause, 
in principle, time-varying phenomena that are observable. 

With a sign reversal of ST a, equation (27) gives a time 
constant ta for exponential growth of /3, 



ta 87rnE(GM)i/2rf/^' 

Evaluation for 4U 0142-1-61, with the parameters of Wang et 
al and a local density E = 300 g cm^^, gives ta ~ 1.1x10^^ s, 
which is too long to be of observable interest. But the value 
of jCw is small for this neutron star and the inner radius of 
the disk is large. Under the assumption that, for example, 
PSR 1828-11 has a disk external to its light cylinder we 
find, given the ATNF parameters for this pulsar, that growth 
times could be quite short, ta ^ 7 x 10^ (n/RLc)^^^ s. 



4 CONCLUSIONS 

We have found that the azimuthal component of the pulsar 
wind momentum-density exerts an r-dependent torque on 
a thin dust or debris disk external to the light cylinder. 
Both the Poynting and relativistic particle components of 
the wind are responsible for ablation of the disk. For disks 
that are optically thick to relativistic protons, the power 
input, which we assume determines the infra-red luminosity, 
is given by equation (15). This result neglects any albedo, 
principally, in this case, reflection of E and B fields at the 
disk surface which could be estimated only by a much more 
complete treatment of this problem than we have attempted. 
The rate of ablation is given by equation (16) in terms of 
the wind luminosity Cw and the cascade parameters £p and 
Af for which we have only qualitative estimates. Under the 
very loose assumption that f2o sin /3 is a universal constant 
for all isolated neutron stars that are formed with dust or 
debris fall-back disks, equation (18) shows that rfS should 
also be a universal constant. From the parameters given for 
4U 0142-1-61 by Wang et al, its value is 1.3 x 10^* g cm. 

It is obvious that our assumption of a common format 
tion spin angular frequency Qo for all neutron stars should 
not be taken too seriously. But if the 4U 0142+61 analysis 
of Wang et al and our analysis of disk ablation are both 
accepted, equation (18) makes it possible to infer what disk 
parameters are possible in other isolated neutron stars and 
pulsars. The conclusion we reach on this basis is that long- 
lived disks of the kind considered here could bo present in 
many pulsars without exceeding published limits on infra- 
red luminosity. 

It has not been possible to answer with any certainty the 



question of whether or not the torque components calculated 
here could be responsible for any observable time-varying 
phenomena. PSR 1828-11 is an interesting case of a pulsar 
that appears to have periodic variations in timing residuals 
and in pulse shape consistent with small-amplitude Eulerian 
precession of its spin a^xis relative to coordinates fixed in 
the star (Stairs, Lyne & Shemar 2000). If the inferred 500d 
period were equated with the growth time ta estimated at 
the end of Section 3.3 for counter- alignment of the inner 
edge of a disk, the required radius would be ri = 2.1Rlc but 
the mean disk density would have to E « 2 x 10^ g cm~^. 
An almost identical value would be required in the radio 
pulsar, PSR B1642-03, which also appears to have periodic 
timing residuals (Shabanova, Lyne & Urama 2001). This 
might be thought a high density for a nominally thin disk 
but, assuming, for example, an outer radius To = 2ri, the 
mass is only Md = 1.5 x 10"'' M©, well below observed upper 
limits. In principle, movement of the ionized inner edge may 
have the potential to produce small changes in the observed 
radio- frequency emission pulse shape and in the pulsar spin- 
down rate. But it is not obvious why the changes should be 
periodic, as observed. 

The more recent case of the isolated neutron star RX 
J0720. 4-3125 in which the inferred blackbody temperature 
and certain other pulse characteristics have been interpreted 
as varying sinusoidally is rather different (Haberl et al 2006; 
but see also van Kerkwijk et al 2007 for an alternative 
analysis). This neutron star with a period of 8.39 s and a 
high magnetic field has some properties in common with 4U 
0142-1-61. The small luminosity, £^ = 4.7 x 10^° erg s"' 
(Manchester et al 2005), is such that, even for n = Rlc, 
the variation-time estimate obtained from equation (28) is 
several orders of magnitude too long to be of interest. 
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